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Particle-particle interactions are of paramount importance in every multi-body system as they 
determine the collective behaviour and coupling strength. Many well-known interactions like electro¬ 
static, van der Waals or screened Coulomb, decay exponentially or with negative powers of the parti¬ 
cle spacing r. Similarly, hydrodynamic interactions between particles undergoing Brownian motion 
decay as 1/r in bulk, and are assumed to decay in small channels. Such interactions are ubiquitous 
in biological and technological systems. Here we confine two particles undergoing Brownian motion 
in narrow, microfluidic channels and study their coupling through hydrodynamic interactions. Our 
experiments show that the hydrodynamic particle-particle interactions are distance-independent in 
these channels. This finding is of fundamental importance for the interpretation of experiments 
where dense mixtures of particles or molecules diffuse through finite length, water-filled channels or 
pore networks. 


Molecular diffusion inside channels and pores is rele¬ 
vant for a wide range of phenomena in natural systems [T] 
as well as technological applications 00 - For example, 
biological channels that transport antibiotics are under 
intense investigation due to their importance for drug 
transport j4]. Diffusing molecules inside protein chan¬ 
nels are closely confined, leading to single-file diffusion 
with hydrodynamic interactions playing a role even at 
these nanometre length-scales M. Transport through 
protein channels can be mimicked by colloidal particles 
and microfluidics chips [7 that confine particles to ID 
diffusion [8|. Such particles undergo random walks in 
one dimension driven by Brownian motion, where the 
corresponding diffusion coefficient critically depends on 
the geometry of the confinement 00. 

Loosely speaking, these Brownian particles receive mo¬ 
mentum impulses from thermal fluctuations of the sol¬ 
vent molecules m, and their resulting motion displaces 
the liquid around it 0HB. This creates a flow field that 
mediates a long-range hydrodynamic interaction between 
the particles [12 . The strength of this interaction is pro¬ 
portional to the flow velocity which decays with distance 
as 1/r, where r is the separation between two uncon¬ 
strained particles in three dimensions (3D) [15] . Intro¬ 
ducing a geometrical confinement modifies the flow field 
that changes the decay rate. For example, the interac¬ 
tion between particles constrained by two parallel plates 
(2D) decays faster at a rate of ~ 1/r 2 [14] . 

In narrow channels the hydrodynamic interactions 
were previously measured to rapidly decay with particle 
separation dSHU]- Indeed, the steady flow induced by 
particle motion in a channel decays exponentially with 
r/2R, where 2 R is the channel width [5j [HJ [18j [19] . 
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Consequently, the particle-particle interaction strength 
is expected to also decay exponentially with their sepa¬ 
ration mm- However, recent theoretical investigation 
suggests that interactions could have a longer spacial ex¬ 
tent and a slower decay rate than previously thought m- 
Furthermore, previous experiments did not capture far 
field hydrodynamics because they used microfluidic chips 
with a groove geometry for mimicking channels 0IIMI7]. 
Sedimentation kept the particles from escaping the mi¬ 
crofluidic groove, but the liquid had no such constraint. 
The lack of controlled experiments in ID confinement 
leaves many unanswered questions about the magnitude 
and spacial extent of hydrodynamic interactions inside 
narrow channels 0H- 



FIG. 1. Trajectories of two particles undergoing Brownian 
motion in closed and open channels. The trajectories are vis¬ 
ibly un-correlated in the closed channel, suggesting the parti¬ 
cles move independently. In the open channel the trajectories 
resemble each other, which leads to strong motion correlation. 


In this paper, we present the first measurement of the 
interactions between two Brownian particles inside a fi¬ 
nite narrow channel, that confines both the particles and 
the liquid. Figure [l] illustrates our experiments, where 
the ‘closed’ channel, shown on the left, has only one end 
connected to a bulk reservoir, while the ‘open’ channel, 
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shown on the right, has both ends connected. Here we 
will demonstrate a fundamental difference between diffu¬ 
sion in open channels and in closed channels mmm- 

The experiments are realised using microfluidic lab-on- 
a-chip devices because they allow for direct imaging of 
particle motion. Inside the chip two large reservoirs are 
separated by a membrane containing closed (figure [2^i) 
and open channels (figure^). All channels are 5 — 17 /im 
long and have semi-elliptical cross-sections with a height 
and width of approximately 800 nm that closely confine 
spherical particles of diameter 2 a = 505 =b 8 nm. The re¬ 
sulting particle to channel size ratio is a/R « 0.6, ensur¬ 
ing particles always move in single-file. Two additional 
large connections are positioned ^200 /im away from the 
narrow channels allowing pressure equalisation between 
the two reservoirs. For the fabrication of the chips, 
we use focused ion beam, photo-lithography and replica 
moulding of polydimethylsiloxane (PDMS) [22] . Cru¬ 
cially, the PDMS chips are then oxygen plasma bonded 
onto a glass slide which provides a bottom wall confine¬ 
ment for the channels. Subsequently, we fill the chip with 
the polystyrene colloidal particles (Polysciences Inc.) dis¬ 
persed in a 5 mM KC1 salt solution that limits electro¬ 
static interactions to a few nanometers. 

An assembled chip is mounted onto an inverted, 
custom-built optical microscope with a high numerical 
aperture oil-immersion objective (xlOO; NA 1.4; UP- 
LSAPO). Using holographic optical tweezers [23h25] we 
position two particles inside the channel, and then turn 
off the laser trapping which releases the particles and al¬ 
lows them to diffuse freely. Their motion is recorded us¬ 
ing a CMOS camera (DMK-31BF03, Imaging Source) at 
a rate of 30 frames per second, until one particle escapes 
the channel. Afterwards, the particle trajectories are ex¬ 
tracted from the images using standard image analysis 
techniques [26]. See supplementary videos SI and S2. 

For data analysis, we divide the trajectories into dis¬ 
placement steps between consecutive frames: Ax\ and 
A ^2 denote the displacements of the first particle (one 
on the left) and the second particle (one on the right), 
respectively. The two-particle interaction strength is ex¬ 
pected to be a function of their separation, and therefore, 
we group the pairs of displacements according to distance 
between centres of the particles. 

For the closed channel, figure shows the displace¬ 
ments of the second particle (Aa^) as a function of dis¬ 
placement of the first particle (Ax\). The distribution 
is circular with points distributed equally in each quad¬ 
rant, suggesting the two particles move independently 
from each other. This is expected for a large particle 
separation of 10a pa ng. In contrast, the distribution 
for the open channel, shown in figure [2ji, is elliptical with 
the major axis along y = x and ellipticity of 0.75. This 
implies that the two particles move in the same direction 
more frequently than in opposite directions, suggesting 
the presence of long-ranged interaction between the par¬ 
ticles. 
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FIG. 2. Comparison between a closed channel (left column) 
and an open channel (right column) experiments. The top 
row (a,b) shows bright field images of microfluidic channels 
containing two colloidal particles. Scale bars indicate 2 pm. 
The middle row (c,d) shows displacements of the second par¬ 
ticle (y axis) as a function of displacement of the first particle 
(x axis). Each distribution contains 1000 displacement pairs 
with an initial particle separation of approximately 2.5 /im. 
The overlaid lines indicate contours for a and 2a from a 2D 
normal-distribution fit. Insets illustrate the direction of mo¬ 
tion for each quadrant. The bottom row (e,f) shows the cor¬ 
relation coefficients versus the separation between the two 
particles. The solid lines indicate the fits to the phenomeno¬ 
logical model. The two-particle motion is strongly correlated 
in the open channel, suggesting a presence of long-ranged non¬ 
decaying interaction. 


We quantify the interaction strength using the Pearson 
pro duct-moment correlation coefficient defined as m ■■ 
p = cov(Aaq, AX 2 ) / a (Ax i) a (Ax 2 ), where cov is the 
covariance, and a is the standard deviation. The value 
p = 0 indicates independent particle motion and p = 1 
corresponds to fully correlated motion. 

Figures^ and[2][ show the correlation coefficients as a 
function of particle separation. Evidently, the correlation 
is stronger and has a longer range in the open channel. A 
detailed examination of the closed channel results shows 
a high correlation when particles are close to one an¬ 
other that exponentially decays to zero in a separation 
of ^4a, as expected. Therefore, we fit a phenomenolog¬ 
ical model, y = A exp (—x/B) + (7, that captures the 
decay rate and has an additional offset parameter. The 
fit yields B = 0.15 ± .02 /im, C = 0.005 db .004. We have 
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introduced the offset to characterise the novel behaviour 
observed in the open channel. As evident from the data, 
the correlation coefficient exhibits the same initial expo¬ 
nential decay, but in stark contrast, it asymptotes to a 
constant offset. This finite correlation coefficient is cap¬ 
tured by the fit to the phenomenological model, yielding 
B = 0.14 =b .02 fi m, C = 0.419 =b .005. This constant, 
non-decaying component is not expected. Furthermore, 
the correlation persists even at the largest measured dis¬ 
tances, suggesting that the two particles interact when 
both are inside the open channel. To the best of our 
knowledge, this is the first observation of such distance- 
independent interactions between Brownian particles. 
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FIG. 3. Correlation coefficients for two interacting particles 
for different channel lengths. The three curves correspond 
to open channels with different lengths, from top to bottom: 
5 /im, 10 /im, 17 /im. The solid lines show the fit to a model 
y — A exp (—x/B)+C, where the offset coefficients were (from 
top to bottom): C = 0.420 ± .006; 0.250 ± .005; 0.101 ± .004. 
Evidently, the long-ranged two-particle correlation coefficient 
decreases with the channel length. Inset illustrates the pro¬ 
posed model for long-ranged hydrodynamic interaction. 

We further investigate the effect of the channel length 
on the particle-particle interaction strength. Figure [3] 
shows the correlation coefficients for open channels of 
lengths L = 5, 10 and 17 /im. The data clearly indi¬ 
cates that the interaction strength decreases with L. In 
the longest channel, shown as circles, we observe parti¬ 
cles interacting at separations of more than 40a. This is 
the largest relative distance measured between interact¬ 
ing Brownian particles, even exceeding the ^20 a separa¬ 
tion measured in bulk Hang. This is a surprising result 
because the geometric confinement typically reduces the 
maximum interaction distance mmm- 

Based on our observations we propose a hypothesis for 
the distance-independent interaction mechanism, that is 
schematically illustrated in the inset of figure [3] Suppose 
the first particle, shown on the left, moves to the right 
due to a thermal momentum impulse. At steady state its 
motion induces a flow that is constrained by the chan¬ 
nel geometry to flow either around the particle or along 
the channel. Importantly, the latter flow has been previ¬ 
ously neglected [sum ng ns] because the studies consid¬ 


ered infinitely long channels, in which the finite pressure 
exerted by the moving particle cannot displace an infi¬ 
nite liquid column. In contrast, we argue that for finite 
open channels flows extend throughout the whole chan¬ 
nel. Consequently, the induced flow along the channel 
has a constant mean flow velocity that is proportional to 
the driving force, i.e. the first particle’s velocity mm- 
This flow encounters the second particle, shown on the 
right, and exerts a viscous drag force on it. The mag¬ 
nitude of this force is a function of flow velocity [30] 
and thus also a function of the first particle’s velocity, 
but is independent of the particle positions. The oppo¬ 
site case of the second particle moving can be inferred 
by symmetry, and the final interaction is a combination 
of the two cases. This gives rise to the particle-particle 
interaction that we observed with our correlation coeffi¬ 
cient measurements and also explains why the value is 
distance-independent in open channels. In a limit, where 
the particle size matches the channel, this problem resem¬ 
bles two pistons in a pipe and the trajectories would be 
perfectly correlated [T2] . On the contrary, in the closed 
channels, the dead end blocks the flow along the channel, 
thus eliminating the non-decaying hydrodynamic interac¬ 
tion. 

We now use this hypothesis to construct an analytical 
model that sheds light into the physics of the interaction. 
The model focuses on physical scalings and omits detailed 
numerical pre-factors from the presentation. Our goal is 
to estimate the typical mean flow velocity in the channel 
resulting from the motion of one particle. 

Consider a spherical particle of diameter 2 a located 
in the centre of an open cylindrical channel of radius 
R and length L. We assume that the particle moves 
instantaneously to the right with velocity U. The fluid 
displaced by the particle must either be pushed along 
the channel to the right channel end, with more fluid 
drawn in at the left, or leak from right to left through the 
thin gap between the particle and the channel wall. The 
pressure increase across the sphere, Ap, is proportional 
to the flow rate in the channel, Q, according to Poiseuille 
law 


where (i is the dynamic viscosity. We neglect any hydro- 
dynamic resistance due to the recirculation from the exit 
of the channel to the entrance, equivalent to imposing a 
periodic boundary condition. Mass conservation around 
the moving sphere (in the frame moving with the sphere) 
leads to 

Q - UR 2 - Rq, (2) 

where q is the leakage flux through the thin gap between 
the sphere and the channel. It is approximately given by 
lubrication theory [31] 
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q ~ —Uh — 


A ph 3 


( 3 ) 


where h = — a is the minimum gap width, and 

k ~ (ah) 1 ' 2 is the characteristic lubrication length-scale. 
Combining these equations and taking the limit h <C R, 
we obtain the flow rate 


Q~UR 2 


f R 3 a V 2 \ 

lw/2 + L/l 5 /V ' 


( 4 ) 


Combining with the mean flow in the channel, 
(u) ~ Q/R 2 , we get the final expression 


M -R 3 q 

U ~ R 3 a + L h 5 / 2 ai/2 1 j 

Our result shows that the mean flow in the channel 
is non-zero, but a function of the channel dimensions R 
and L. Importantly, our model reproduces the observed 
decrease with the channel length - scaling with 1/L. The 
long channel limit (L oo) gives no external flows (Q = 
0) equivalent to the closed channel case. 

In order to quantitatively capture the experimental re¬ 
sults, we employ numerical simulations. Using finite ele¬ 
ment analysis software (COMSOL Multiphysics v4.4) [Qj 
we solve the Stokes equations for two spherical parti¬ 
cles inside a cylindrical channel with a diameter chosen 
such that the areas match the experiment parameters. 
Periodic boundary conditions are set on the ends while 
all other surfaces were set to no-slip boundary condi¬ 
tions. We apply an instantaneous velocity on one par¬ 
ticle and calculate the resulting drag forces on both par¬ 
ticles. A velocity independent measure of these forces 
is known as a friction matrix (£) [29l p. 226]. For 
the two particle system in ID it is a 2x2 symmet¬ 
ric matrix, where the diagonal terms (Ci,i 5 (2,2) de¬ 
scribe hydrodynamic drag experienced by the moving 
particles (first; second), and the off-diagonal elements 
(Ci,2 = C2,i) correspond to a force exerted by the mov¬ 
ing particle on the other particle. The friction matrix 
values are computed directly from the numerical sim¬ 
ulation. In addition, the friction matrix can be ob¬ 
tained from the experiments using a diffusivity matrix, 
D = /c#TC -1 [29], with components D i:i = (Ax 2 )/2At 
and Dip = D 2 ,1 = (Axi Ax 2 )/2A t [32] . This allows us 
to compare the hydrodynamic interactions predicted by 
our simulation with the measured data. 

Figures [4^i and illustrate the typical flow patterns 
computed with our simulation. Notice that the flow in 
the closed channel curls around the particle and does not 
extend far into the channel. This contrasts with the open 
channel, where the flow around the particle is weaker and 
there is a Poiseuille flow along the whole channel length. 
This moving fluid column exerts a force on the second 


particle that gives rise to the non-decaying interaction. 
Figures H and [4]l show a quantitative comparison be¬ 
tween the numerical simulation and experimental data, 
where the solid lines indicate the simulation results and 
the points indicate the experimental values. The fric¬ 
tion with the channel walls is the same for both parti¬ 
cles, leading to the overlapping curves ((j ? i and £ 2j2 ). 
Meanwhile, the interaction force (£ 1 , 2 ) asymptotes to a 
non-zero value for the open channel only, similar to the 
correlation coefficient that was reported above. At small 
separations the discrepancy between simulation and ex¬ 
periment is likely caused by electrostatics and finite frame 
rate. Meanwhile, at large separations the values compare 
very well, with the largest discrepancy below 20%. This 
is very good agreement given the approximations made 
about channel shape and width. 

The agreement between the numerical simulation and 
the experimental data suggests that the proposed hy¬ 
pothesis captures the physics of non-decaying interac¬ 
tions between the Brownian particles. Crucially, this in¬ 
teraction requires flow through the ends of the channel, 
that is attained by keeping the two reservoirs at equal 
pressures. In our experiments, the pressure equalises via 
liquid recirculation through the secondary large chan¬ 
nels placed far away from the narrow channel. In bi¬ 
ological cells, this can happen through the membrane 
or through dedicated water transport protein channels 
known as aquaporins pQ. In porous materials, the pres¬ 
sure can equalise through other interlinking channels [2]. 

Our observations have far-reaching implications for dif¬ 
fusion processes inside channels. One of the predictions 
from our model is that the viscous drag on particles is 
smaller than the theory expects for narrow channels (fig¬ 
ure S5). We therefore performed an additional experi¬ 
ment that measured the diffusion coefficient of a single 
particle inside an L = 10 pm open channel (method [9]). 
The measured diffusion coefficient is indeed 40% higher 
than expected [5j |TT] : D x /Dq = 0.126 =b .006, where 
Do is the diffusion coefficient in the bulk (figure S6). 
This potentially explains the experimental discrepancies 
noted in previous studies for channels with similar di¬ 
mensions [SESIEI]- Because the diffusion coefficient is 
important for predicting transport rates, our results di¬ 
rectly impact channel transport models [35]. Moreover, 
the long-ranged interaction between particles should lead 
to a cooperative behaviour that could enhance transport 
across channels. Also, since the interactions are not lim¬ 
ited to two particles, they should persist for channels 
filled with three or even more particles. This makes 
our results relevant for polymers in confinement because 
the monomers can interact non-locally [36], and also for 
single-file systems because the particles interact not only 
with the closest neighbours [8]. 

In conclusion, we utilised a highly controlled microflu¬ 
idic system coupled with holographic optical tweezers to 
investigate the interaction of particles in confinement. 
Our measurements prove that interactions extend over 
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FIG. 4. A comparison between experimental data and simulation results. Top row illustrates flows inside the (a) closed 
and (b) open channel based on simulation results. The proportions were chosen for illustrative purposes. The bottom row 
compares the friction matrix values as function of particle separation for the (c) closed and (d) open channel. The points 
indicate experimentally measured values, while solid lines show the corresponding simulation results obtained without any 
fitting parameters. All values are normalised with the friction coefficient in the bulk: £o = 67r fia. Error bars are shown on every 
fifth data point. Arrows highlight the non-zero interaction term which agrees well with our simulation results, suggesting that 
the non-decaying interaction is caused by the induced flow inside the channel. 


the full channel length and have a constant strength 
that does not decay with particle separation. We ex¬ 
plain the coupling mechanism using hydrodynamics with 
both a qualitative analytical model and quantitative 
comparisons with the numerical simulations. The ex¬ 
cellent agreement between the theory and experiments 
suggests that we fully understand the properties of hy¬ 
drodynamic particle interactions in microfluidic channels. 
The non-decaying interaction extending throughout the 
whole channel has important implications for the mod¬ 
elling of transport through channels as well as for the 
interpretation of experiments investigating particles dif¬ 
fusing in close confinement. 
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